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Abstract 

The bilocal effective action is developed as a viable tool to study 
mesons in the large-N limit of QCD. Several results from current alge- 
bra are derived in a new way using this action. The question of gauge 
invar iance is discussed in an appendix 



1 Introduction 



Although quantum chromodynamics has been extremely successful at de- 
scribing quarks and gluons perturbatively, it has yet to make many firm 
predictions about mesons and baryons. How can the composite states be de- 
scribed, starting from the Lagrangian for the constituents? Many attempts 
invoke the constituent quark model and/or the potential or bag model [0], 
none of which can be derived from QCD. The Nambu-Jona-Lasinio model 
offers another approach at the expense of approximating the gluon prop- 
agator by a delta-function. 

Perhaps the most puzzling feature of these models is that they work so 
well. Surely these models' successes, in spite of gross dynamical approxima- 
tions, indicate that each model shares with QCD a common symmetry (or 
other fundamental principle) that is actually dictating the favorable results. 
Those predictions that depend explicitly on the dynamics are (perhaps) the 
less successful ones. However, it is impossible to tell which of the mod- 
els' predictions are artifacts of the dynamical assumptions, and which result 
from the fundamental principles, since the dynamical approximations are 
employed from the beginning. 

It therefore seems worthwhile to attempt a systematic approach to de- 
riving hadron properties from QCD, making as few assumptions as possible. 
To do this, I have concentrated on the mesons, and have used the effective 
action for bilocal operators. 

However, some approximation must still be used. I have used the large 
number of colors (large-N) expansion, which offers a formal justification for 
the "quenched" approximation. In this scheme, a meson consists of a single 
quark-ant i- quark pair (and an indefinite number of gluons), other virtual 
quark pairs being suppressed by powers of 1/N [|J. Therefore the meson n- 
point functions contain a single fermion loop; this simplification allows one to 
write down — at least formally — exact expressions for these functions. And 
even though the n-point functions still cannot be computed, it is possible to 
derive relations between them. 

In this paper, the effective action for large-N QCD is derived. The vacuum 
is discussed first, then the mesons. The approach illuminates the dual nature 
of the pion: it is both a bound state solution to the Bethe-Salpeter equation, 
and it is a Goldstone boson, present because the vacuum state spontaneously 
breaks chiral symmetry. Finally, heavier mesons are considered. A little- 
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known sum rule is derived in a new way using this approach; the sum rule 
agrees quite well with experiment. The appendix addresses the question of 
the gauge invariance of the results. 



2 The Effective Action 

We begin with the formalism of Cornwall, Jackiw, and Tomboulis ||. They 
derive an effective action for the bilocal operator G a b(x, y) = {Tip a {x)'4)i > {y)) ■ 
The indices a and b give the spin, isospin, and color of the fermion; they are 
suppressed below. The effective action is 

r[G] =itr log[G] -tr[(i$-m)G]+r 2 [G] (1) 

The matrix m gives the current quark masses (~ 5 to 10 MeV for the u and 
d quarks). r 2 [G] is the sum of the two-particle irreducible (2PI) diagrams, 
such as in Fig. 1. That is, it is the sum over all diagrams that remain 
connected after breaking any two fermion lines. The fermion propagator for 
these graphs is equal to G. T 2 can be expressed formally as 

•2 -3 

iv 2 = l -{{T$4i>){Um - ^{{Tm^m^mm + ■■■ 

where the double brackets mean that the path integral over all possible field 
configurations must be performed, and in the end only the 2PI diagrams are 
retained. The sum over possible A^'s is subject to a gauge-fixing constraint; 
the appendix shows that physical observables are, however, gauge invariant. 

To leading order in the 1/N approximation, there is only one fermion 
loop in each diagram 0]; therefore, replacing all fermion propagators by 
G, including the appropriate symmetry factors, and remembering an overall 
minus sign for the fermion loop, 

i 2 i 3 



iT 2 [G} = 

= <(trlog(l + zG4)» 



2 tT((G4GA)) - -tr((G4G4G4)) + 



((tr logzlG--l» + ((tr log4» 
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Therefore, ignoring any constants independent of G, the total action is 



i 



T[G] = itr log[G] -tr[(i^-m)G] -i((tr log ^G - -j 1» (2) 

3 The Ground State 

The function G satisfies the equation of motion 5T/5G = 0: 

jG?- 1 -(^-m)-i« ?7 i- r )> = (3) 

^ _ J 

The rest of this section does no more than define some new notation so that 
Eq.(|3|) can be written more conveniently as Eq. fl7|), below. 

In the ground state, translation invariance requires that G(x, y) be a 
function of x — y only: G(x,y) = G (x — y). Using parity conservation 
and Lorentz invariance, Go can be expressed in terms of p, the momentum 
conjugate to x — y, as 

•Go 1 = <P 2 )t ~ KP 2 ) 

for two unknown functions a and b. If b ^ when the quark mass m = 
0, the propagator Go spontaneously breaks chiral symmetry. Coleman and 
Witten have shown || that large-N QCD does in fact break chiral symmetry 
spontaneously, so we expect that b ^ 0. 

Defining a = a/ (a 2 p 2 — b 2 ) and /3 = bj (a 2 p 2 — b 2 ), we can write 

G = i(a$ + P). 

Then Eq. (|3]) becomes 

(a-l#-(6-m) = (( _\ )) (4) 

a r ^ + P 

The notation can be further condensed by defining f = aft — 1/jjL. Then 
using the identity 

'Tf/^r 1 ) (5) 
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and equating an even number of gamma-matrices on both the left- and right- 
hand sides of Eq. (||) , we get 

b + m = -(( „_ 1 r(3rip ^r 1 )) (6) 



This expression is simplified by using Eq. ([5]) to write 
1 
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Then Eq. (|6]) becomes 



Because this equation gives a consistency condition for the order parameter 
b that appears as a result of the spontaneous breakdown of a symmetry, it is 
analogous to the gap equation in the BCS theory of superconductivity. 



4 Excited States 

The previous section concerned the vacuum state where G(x, y) = Go(x-y). 
Let us now look for the excited states of the theory. They can be found by 
allowing G(x,y) = Gq(x — y) + 4>{x,y) and requiring that the total G still 
satisfy the equation of motion, ST/SG = 0. To see that this procedure is 
correct, note that it implies that 



5T 



6T_ 

6G 



+ 



5 2 T 



Go 



5G 2 



Go 



The first term on the right hand side is zero by definition; using Eq.(|I|) the 
second term gives 

c2 i 

2 



— iGq ^c/)Gq ^ 



5 2 r 2 



5G 2 



0. 



(8) 



G„ 
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Since 1^ is the sum of 2PI graphs, 6^2/ SG 2 is the irreducible four-point 
function. Therefore Eq. (J8|) is simply the Bethe-Salpeter equation for 0. 
This confirms the suspicion that represents an excited state: it is the 
Bethe-Salpeter wavefunction [0] of the quark-anti-quark bound state. 

We can now write the effective action for the meson fields by expanding 
T[G] around G : 

T[0] = -~tr U^-^r) + ~tr U^^^r 

+ ^((0-^0—^)) - l t r<<0— ^0— ^0— L^) + . . . (9) 
The equation of motion for is 

TT^TT = (( r ! i $ r \ i )) + (interactions) (10) 

5 Pions 

Armed with this general formalism for mesons, we can look for particular 



solutions to Eq. fllPl) . It is easiest to start with pions, since we expect 
them to exist on general grounds due to the spontaneous breakdown of chiral 
symmetry. With some foresight, then, we make the ansatz 

<P = X {x-y)W{y) 1 ,T a (11) 

for the pion fieldQ (T a is a generator of the flavor group.) We expect that 
H a (y) oc exp(iky), where k^ is the four-momentum of the pion j7|. For 
simplicity, consider the case when k^ — > (which is only possible for physical 
pions if the mass of the pion is zero). Then Eq. (|T0|) becomes 

{aj> - b)xi,T a {ai> -b) = ((^L_ xl5 T a ^)) 

Anti-commuting the 75 through, and using the fact that x commutes with 
p M , gives 

-y-_ i2 >* = «7*b*^» (12) 

1 As Bethe and Salpeter point out 0, we could use 4> = xi x ~ y)Tl a {ax + (1 — a)y)j5T a 
for any a between and 1. The value a = is chosen simply for convenience later. 
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If x — ft, this equation is exactly satisfied due to the gap equation, Eq. (|7|), 
when the quark mass is zero! The fact that a solution is possible when — > 
means that indeed the pion is massless.Q And this approach clearly shows the 
dual nature of the pion as both a bound state and as a Goldstone boson: Eq. 



(|12D is the Bethe-Salpeter equation for a meson; yet it is identical to the gap 
equation for the vacuum,^ Eq. (^). Furthermore, the internal wavefunction 
of the pion can be related to a matrix element of the vacuum (for rif flavors 
and N colors; tr' is a trace over all indices except spacetime): 

xO -y) = -y) = -; — T7 tr '( G o) 

4rifi\ 
tr'(0|V(aM2/)|0>. 



An f N 



The effective action for pions can now be constructed, using the ansatz of 
Eq. (|TT|) with x — P, and relaxing the condition H a (y) = const. The quark 
masses are assumed to be nonzero; for simplicity only the two-flavor case, 
with m u = md = m i will be considered. Defining II = IT a T a , 

r = \f (Ptji(at-b)prin(at-b)) 

-^ mb n 7 L- p ^u f L^)) + .... (is) 

The action above can be used to determine the mass of the pions in the 
presence of explicit chiral symmetry breaking (i.e. when the quarks have a 



small mass). One can perform a derivative expansion on the field II |T1 
commute the 75 matrices through other quantities until they stand next to 
each other, and use the definition of j3 in terms of a and b to obtain: 

r = -^tr(&/m 2 ) + ^<<^/?^ym 2 >> + ^( W 2 + . . . 



2 It would appear that the same argument should hold for the U(l) meson, the 77'. 
However, for this meson the 1 /N corrections connected with the anomaly cannot be ignored 
§■ 

3 The equivalence between the Bethe-Salpeter equation and the gap equation was also 
shown, for the case of scalar exchange between fermions, by Migdal and Polyakov [ fio| . 
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where the constant Z can be determined from the derivative expansion of 
Eq. flHf). The gap equation, Eq. (0), further simplifies the action: 

1 1 

T = — tr(m/m 2 ) + — {d^Wf + (interactions) 

2 2Z 

The physical pion field, 7c p h ys , should have the canonical normalization for 
its kinetic energy. Therefore the field IT is related to 7r p h ys by Tt p h ys = 11/ vZ. 
Then the mass of the pion can be read off from the action above (using 
tr (T a T b ) = \5 ab ): 

2 imtr(S)Z 
m « = 4 • 

The normalization constant Z can be evaluated in the following way: 
the Bethe-Salpeter amplitude is = (0\ip(x)ip(y)\ii a ) = /3(x — y)Il(?/)75T a . 
Therefore 

(0|V(0) 75 TV(0)|vr a ) = -^ afe 4iV/3(0)n(0) 

= --8 a hr((3)\fZ 
On the other hand, current algebra can be used to show that 
(0|^(0) 75 TV(0)K> = ^-{^)5 ab = ^-5 ab . 

Therefore \[Z = — 2//^, and 

2 .. m($ifj) 



f 2 



The effective action therefore reproduces this classic result for the mass of 
the pion ||. 



6 Other Mesons 

Encouraged by the effective action's ability to reproduce pion results, it is 
time to move on to heavier mesons. It is possible, with one additional as- 
sumption, to produce another current algebra result relating the mass of the 
p to the mass of the a%. 
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The p and a\ mesons can be parametrized as 

p = ({x- y)R;(y)rT a ; &1 = ((x - y)A«(y)n 5 T a 

The extra assumption is that the internal wavefunctions of the p and a± are 
nearly equal: £ ~ £. The relation between the p and ai can then be sketched 
as follows: using the same method as for the pions, the part of the effective 
action that determines the mass of the p contains a term 

where Z p is the normalization constant for the p-meson. Similarly, the mass 
of the ai contains a piece 

1 2 i 1 . 1 



-m 



tr((^ 75 — — £ 7m 7 5 — — )) 



2Z a ai 32 '°f+P f+/? 



Therefore, 



z a z p ie trU57 r+pr 1,i f+p" 

11 1 

= ^^T^ 75 ^^ 75 ^ (14) 

Notice that the last term has the form of a p — it — a>i interaction (cf. Eq. 
(H)). Using the complete expressions for the masses of the p and a\ confirms 
this argument. Defining the coupling f apn by the interaction 

we find 

9 2 f 

Wilt ™n Tnnir 



Z a Z p JZ n Z n Z« 
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This sum rule can be put into a more convenient form by using the following 
definitions of g a and g p , 



(o\miuT b m\p;) = g P ^g^ = - Vvmo jz, 
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(0|^(0)7,75TV(0)|< M ) = g a 5 ab g pu = -±5 ab g,M0^fz~a 
and using \fZ~ K = —2/f n to obtain 

m ai _ m p _|_ f-wfapir (15) 

gl g 2 P g a g P 

This sum rule was originally derived by Schnitzer and Weinberg |I5| using 
current algebra. In order to test the relation numerically, the masses of the 
p and a\ can be taken from the Particle Data Group [f[^| ; f n « 92 MeV 
g p can be extracted from p — ► e + e~ decay, g p = 0.117 ± .005 GeV 2 0; and 
/ ap7r can be obtained from a\ — > pn decay, f ap7r = 4.6 ± .2 GeV [|y|. The 
only problematic quantity is g a . Assuming that the decay t —> pn always 
proceeds through the ai resonance, g a = 0.177 ±0.014GeV 2 pi. However, 
Weinberg's sum rule [ |H}] , 



g^ , 

2 T~Kt 2 



gives a value g> = 0.15GeV 2 . Using the first value of g a , Eq.(|l5|) gives 

48GeV~ 2 = 64GeV" 2 , 
which is about a 25% discrepancy. However, the second value for g a gives 

67.2GeV" 2 = 67.4GeV" 2 



7 Conclusion 

The utility of the effective action is that it provides — as do all effective 
actions — a convenient bookkeeping device for the quantum numbers of the 
particles. However, unlike the purely phenomenological actions, the terms 
in this action automatically appear with the correct symmetry properties, 
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since they are derived formally from QCD. Furthermore, the bilocal action 
provides a systematic approach to the study of mesons, housing the gap 
equation, the Bethe-Salpeter equation, and scattering amplitudes all under 
one roof. It also provides a starting point for more detailed model-based 
calculations — if one would like to use a particular ansatz for the Green func- 
tion Go, for example. However, as the above results demonstrate, a number 
of calculations can be performed even in the absence of such assumptions. 
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Appendix: Gauge Invariance 

This appendix addresses the difficult question of gauge invariance. It turns 
out that the effective action is explicitly gauge dependent, as is the Green 
function Gq. However, it will be shown that meson masses, as well as physical 
scattering amplitudes, do not depend on the gauge. The proof depends on 
a generalization of methods used for the scalar effective potential ||16|| . An 
excellent pedagogical introduction to the subject is given by Ref. fT 



We begin with the QCD path integral, with various source terms added. 
We first add a bilocal source K\ it will be replaced by G in the Legendre 
transform. For convenience, we also need a source h that is independent of 
x, and sources Q and R, which are bilocal. To condense the notation, the 
spacetime dependence of the variables, when expressed, will be written as 
indices rather than as arguments (e.g. Q xy instead of Q(x, y)). The path 
integral, including the gauge-fixing term and the ghosts, is 

Z[K,h,Q,R] = J DiljDijjDriDfjDA^expi^J dx - 4 - m)ip - -^F 2 

+ J dxdy {-MQ xy tVyT a - t V a x T a Q xy )^y - ^ x ri a x T a R X y^f y T h ^y - ^ x K xy ijj y ]\ 
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The trick is to use the Ward identity that results from the fact that the path 
integral is invariant under the BRS transformation; this can be converted 
into the statement that physical observables do not depend on the gauge 
parameter a. Under the BRS transformation || (using the ghost equations 
of motion to simplify the result), 

5Z = J D[fx] Q-h(d^) 2 - $[K, irfT a \i\) - 2^{rj a T a Qr] b T b )^j e iS = 

In this equation the variables are treated as having spacetime indices, and all 
indices are summed over. Now the sources have been carefully introduced in 
just such a manner that the above equation can be rewritten using derivatives 
with respect to the sources: 

(^-L +t ^ K ^ )+2tr{Q "ak ) ) z=0 (A1) 

From here the generator of connected Green functions W is introduced in the 
usual way: Z = exp(iW). The effective action T is the Legendre transform 
of W: 

T[G] = W[K] - tr{GK) 

with 

dW dV 
dK~ G ' dG~~ K 
The sources h, Q, and R are not Legendre-transformed. Eq. (|A.l|) becomes 

Eq. O is the main Ward identity from which everything else in this 
appendix is derived. 

To begin, we take the partial derivative of Eq. ( |A.2|) with respect to h, 
set Q = h = 0, and omit terms which vanish due to ghost conservation: 

fdT d 2 T \ dT , A . 

-'idGdhdQ)^^- (A ' 3) 

This equation states that the action changes when a is changed. However, 
the value of T at the point where dT/dG = does not in fact depend on a. 
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To learn about the meson two-point function, we take the derivative 
d 2 1 ' dG xy dG zw of Eq. (|A.3| ). (When G and Q appear below without (space- 
time) indices, the indices of G are contracted with the indices of Q.) 

( d 2 T d d\ d 2 T <9 3 r d 2 T 

+ 01- 



dhdQ dG da ) dG xy dG zw dhdQdG xy dGdG zw 
d 2 T 3 3 T or <9 4 r 

= (A.4) 



8G xy dG dhdQdG zw dG dhdQdG xy dG z 

Two observations simplify this equation somewhat. First, this equation 
should be evaluated at G — Go, where dT/dG = 0, so the last term of 
Eq. ( A.4j ) vanishes. Second, it can be shownP] that the function d 2 T /dhdQ, 



evaluated at G — G , is equal to —a dG /da. Therefore the expression in 
the first term of Eq. (|A.4|) , 



8 d 2 T 8 

a- 



da dhdQdG 



G=G 



( d dG d 
a\ 1 

V da da dG 




gives the total a-dependence, both implicit and explicit. 

Eq. (|A.4 ) shows how the inverse two-point function changes as a changes; 



we also need to know how the Bethe-Salpeter wavefunction changes. We can 
multiply Eq. ( |A.4|) by the Bethe-Salpeter wavefunction (p zw from the right, 
and use the Bethe-Salpeter equation (d 2 T /dG uv dG zw )(p zw = (sum over 
repeated indices!) to get 

/ D d 2 T \ d 2 T d 3 T 

[ a Da 8G xy dG zw ) 02W 8G xy dG dhdQdG zw ^ ~ 

Under a change of gauge, a — > a + 5a, the wavefunction (ft must change by 
a certain amount 5(f) = (D<j)/ Da)5a in order for it to still satisfy the Bethe- 
Salpeter equation. The amount by which it must change can be read off from 
the above equation: 

D d 3 T 

Da dhdQ zw dG z i w i 



The analogous case holds for the scalar potential |T(| [l7|. In this context it can be 



proved by taking the derivative of Eq. (A.l) with respect to K, Legendre transforming 
the result, and setting G = Gq. 
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This ensures that 



a 



D 
~Da 

i> + > 



d 2 T 



OG xy dG 

z 







is the correct Bethe-Salpeter wavefunction 



which means that indeed 
in the new gauge. 

It proves more convenient to express Eq. (|A.5 ) in momentum space. 
Recall that the Bethe-Salpeter wavefunction can be factored: <p xy = xi x ~ 
y)H(y) where \ describes the internal motion, and H the center-of-mass 
motion. Fourier transforming this relation gives <j)(q; P) = x{q)H(P), where 
P obeys the relation P 2 = M 2 . Now writing the Fourier transform of the 
function appearing in Eq. flA.5 ) as 



d 3 T 



dhdQ xy dG zw 
Eq. (|A.5| ) becomes 
D 



dq dk dP' f (3) (q, k- p') e ii^-y) e ik.(z-w) e iP'.(y-w) 



a—( X (q)H(P)) = -J dkT^(q,k;P) X (k)H(P) 

(Here translation invariance has forced P' to equal P.) The function H does 
not change under the gauge transformation, and we have 



Da 



dk T {3 \q,k;P) 



P 2 =M 2 



X(k) 



(A.6) 



At last we are ready to show that meson masses are gauge independent. 
H is a local function, and we can regard T as a functional of H, so that the 
derivatives of T with respect to H give the one-particle irreducible amplitudes 
for meson scattering. To do this, we write H = H(P), where P can now take 
on any value. Then the Fourier transform of the inverse two-point function 
is 

d 2 T „ r d 2 T 

dq dp x (p) (p,q;P)x(q) 



dH 2 



(P) 



dG 2 



Using Eq. ( |A1D and Eq. ( [OP we find that 



d 2 T 



a 



-( 

Da \ dH 2 



P 2 =M 2 , 
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That is, if P 2 = M 2 , the value of the inverse two-point function does not 
depend on the gauge parameter a. Therefore the condition P 2 = M 2 , which 
causes the inverse two-point function to vanish, causes it to vanish in all 
gauges. This means that the value of the (renormalized) physical mass M, 
defined as the point where the inverse two-point function vanishes, does not 
depend on the choice of gauge. 

The three-point functions can be analyzed in a similar way. We take the 
derivative d/dG uv of Eq. ( |A.4|) ; multiply by the wavefunction xiM) where k 
is the momentum conjugate to {u — v); then use Eq. ( |A.6|) to find 



D ( d 3 T \ 

if all three particles are on shell. That is, the off-shell amplitude in general 
depends on the gauge, but the physical three-point scattering amplitude is 
gauge independent. 

The three-point function calculated in section |6] was not taken on shell, 
but with the external momenta going to zero. However, Eq. ([14]) is gauge 
co-variant, by construction. It relates the masses to the wavefunction renor- 
malization constants and the three-point function, all evaluated at p = 0. 
The subsequent derivation of Eq.(^) relies on the additional assumption, 
common in meson physics, that the Z's and f apn depend only weakly on the 
momentum, no matter what the gauge. 

Finally, there is an amusing property of the four-point function, a Ward 
identity for which can be obtained by taking the derivative d 2 /dG uv dG st of 
Eq. ( |A.4| ). Curiously, the one-particle irreducible four-point function is not 
gauge-invariant, even when the external particles are on-shell. Only the com- 
plete four-point amplitude — the sum of irreducible and reducible diagrams — 
is gauge invariant. The fact that the irreducible four-point function mixes 
with the reducible functions under a gauge transformation means that a me- 
son cannot be thought of as a distinct particle that maintains its identity 
throughout the reaction! 
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Figure Captions 



Fig. 1 The two-particle irreducible diagrams contributing to the effective 
action. 
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